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Saturated stellar images on digitised photographic plates are many 
times greater in area than the 'seeing disk' seen in unsaturated CCD 
images. Indeed the flux profile of a bright star can be traced out 
for several degrees from the star's centre ' 2,s . The radius of the sat- 
urated stellar image can often be directly related to the magnitude 
of the star, a fact well known and exploited in iris photometry 4 ' 5 . 
In this work we compare the radial flux profile of stars in the ap- 
proximate range B~9 to ~13 mag, obtained from scans of plates 
from the Bamberg Sky Patrol archive, with a profile of the form 
measured by King 1 . We show that simple aperture photometry of 
saturated stellar images, obtained from photo-positives of scanned 
photographic plates, yield data that are in agreement with simula- 
tions using a (saturated) synthetic stellar radius profile. Raw plate 
magnitudes from this aperture photometry can be easily and satis- 
factorily transformed to standard magnitudes, as demonstrated in 
a recent study carried out by the current authors 6 . 

Introduction 

The radial profile of a star observed by a ground-based imaging system is 
known to be far more extended than would be expected from the instantaneous 
seeing disk or from diffraction effects in the instrument. King 1 , Kormendy 2 
and Racine 3 provided examples of the stellar profile, and discussed possible 
causes. It appears that both instrumental and atmospheric effects are likely to 
be contributing factors 3 . Fig. 1 shows the profile derived by King (open circles) 
for a zeroth-magnitude star from photographic data, together with a polynomial 
fit to the data (solid line) which we will use below as our synthetic profile. We 
have adopted this, rather than the more recent, CCD-derived, profile of Racine 3 
as we are concerned with aperture photometry of scanned photographic plates. 

Significantly, for radii greater than about 10 arc seconds, there is an ap- 
proximate power law relation between flux (in magnitude per arc second) and 
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Figure 1: Stellar profile, in B magnitudes per square arc second, from King 1 . 
Open circles represent data points digitised from King's Figure 1. The line 
represents a polynomial fit to the data (from 0.5 arc second radius outwards) 
used to create the synthetic stellar profiles employed in the current work. 
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radius 7 . Hence, for iris photometry of a saturated stellar image of radius of 
10 arc sec or more, when the star is centred in the iris diaphragm, there is 
a linear relationship between the logarithm of the measured transmission and 
magnitude 5 . There is a direct, but inverse, relationship between the darkening 
of the original photo-negative plate and the brightening of a photo-positive im- 
age. Consequently we might expect, for a digitised photo-positive plate scan, a 
near-linear relationship between the plate magnitude (logarithm of the 'counts' 
from aperture photometry) and stellar magnitudes in this case. 

Aperture photometry of photo-positive plate scans 

We explored the above suggestion using a series of scans made of Bamberg 
Sky Patrol plates (e.g. see Strohmeier and Mauder 8 ), from the plate archive 
of the Dr. Remeis Observatory, Bamberg. The scans were originally made for 
a study of the active star CF Octantis. A short description of the Sky Patrol 
programme and more details of the CF Octantis results are presented in Innis et 
al. 6 . In brief, the Sky Patrol plates were obtained with 0.1-m diameter cameras 
at three Southern Hemisphere sites between 1962 and 1976, and have a typical 
limiting photographic magnitude near 14. An Epson flat-bed scanner was used 
to digitise 375 plates containing the field of CF Oct, at a resolution of 5.25 arc 
sec per pixel. The scans were converted to photo-positives and saved as FITS 
files. The analysis package IRAF was used to obtain aperture photometry of 9 
field stars (ranging from B^9 to 10.8) and CF Oct on 353 plate scans, with the 
other 22 plates being unusable due to cloud, plate scratches, or other effects. 

Example raw plate magnitudes for two field stars, HD 195460 (B^8.99) 
and CP -80 980 (B-10.62) are shown in the top panel of Fig. 2. Both stars 
showed saturated cores on the plates. We used a fixed aperture of 20 pixel 
radius (~105 arc seconds) for the measurements (see Innis et al. 6 for more 
details). Large and sudden vertical shifts in the mean levels correspond to 
annual breaks in the data runs. Data from three different sites are shown. While 
the magnitudes vary significantly from plate to plate, by up to ~2 magnitudes, 
they do so in a reasonably consistent manner. The magnitude differences (lower 
panel, same vertical scale) are much more uniform. For these two stars the 
standard deviation of the differences was 0.16 mag. As noted, 9 field stars (and 
CF Oct) were measured on each plate. For clarity of presentation results for 
only two field stars are shown, but these illustrate typical behaviour. As part 
of a further study, 10 additional field stars (down to B^13) were measured on 
a subset of 90 plates. Similar results were found, albeit with slightly increased 
observational errors. Combining the two field star data sets for the 90 plates 
gives plate magnitudes for 19 stars ranging from B~9 to ~13, which forms the 
basis of our present study. 
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Figure 2: Example sequence of raw, untransformed, plate magnitudes, from 
IRAF aperture photometry, of two stars, measured from ^350 digitised Bamberg 
sky patrol plates from 1965 to 1969. Both stars showed saturated cores on the 
plates. Top panel: raw plate magnitudes of HD 195460 (circles) and CP —80 980 
(+ symbols). The large jumps near sequence numbers 75 and 300 correspond to 
annual breaks in the plate runs. Lower panel: Differential magnitudes between 
these two stars, plotted on the same vertical scale as the upper panel. The 
standard deviation of the differences was 0.16 mag. We measured up to 19 
stars on each plate, but for clarity results for two stars only are plotted here. 
Although the zero points change significantly from plate to plate, by nearly 2 
magnitudes, the magnitude differences remained much more consistent. 
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Synthetic stellar profiles 

Fig. 3 shows a radial plot of an observed ('sky subtracted') stellar profile 
(square symbols), obtained from a digitised plate, plotted as arbitrary intensity 
versus 5.25 arc second pixels. The observed profile is saturated out to about 6 or 
7 pixels in radius. Also plotted is part of the polynomial fit to the King profile 
(from Fig. 1), with a saturation limit imposed so that the radius of saturation is 
6 pixels to approximately match the observed profile, and then scaled vertically 
so that the saturation intensities match. The agreement is reasonable although 
the observed profile falls more rapidly than the synthetic profile at larger radii. 
The difference at large radius is, at least in part, probably a consequence of 
our polynomial fit (sec Fig. 1), which appears slightly high near 100 arc sec 
(~20 pixels). For the purposes of this study we believe the small difference can 
be neglected, as we are attempting to illustrate a general principle rather than 
derive an exact match to the stellar profiles recorded on the Bamberg plates. 

We generated a series of synthetic stellar radial-flux profiles, based on our 
fit to the data of King 1 , spanning a range of 4 magnitudes. These are shown in 
the left panel of Fig. 4 on a semi-log plot, with the x-axis in arc seconds. We 
then applied a common 'saturation limit' to all profiles to mimic the saturation 
level of the photographic emulsion, and resampled the data to a 5.25 arc second 
pixel size. The right hand panel of Fig. 4 shows the resulting profiles. The 
radius of saturation varies from around 2 to 13 pixels. As explained below, the 
saturation intensity level was chosen so that the resulting radii of saturation 
matched the observed values from the plate scans. We stress that the only 
adjustable parameter in this matching is the choice of saturation intensity level, 
and this is applied to all the synthetic profiles. 

Radii of saturation 

Fig. 5 shows, as open circles, the measured radii of saturation for our 19 
field stars, ranging from B^9 to ^13. We measured radii on 6 digitised plates 
whose plate magnitudes were close to the ensemble mean determined from 90 
plates, and then took averages for each star. This procedure should minimise 
possible changes in stellar profiles caused by instrumental changes made during 
the survey. The calculated radii of saturation of the synthetic profiles of Fig. 4 
are represented by the solid line. We adjusted the common saturation level for 
the synthetic profiles to get a reasonable match (by eye) between the calculated 
and measured radii of saturation. As noted, this is the only parameter that we 
adjusted, after fixing a zero-point for the synthetic magnitudes. The agreement 
is reasonable for the larger radii, but is poor for the fainter stars. In part 
this may be a consequence of the limited resolution of the plate scans: The 
synthetic profiles were calculated on a 1-arc second axis, then resampled to a 
5.25 arc second pixel scale to match the plate scan resolution. When the area 
of saturation of the profile on the plate is comparable to the 5.25x5.25 arc sec 
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Figure 3: Radial profile plot showing a comparison between an observed (sky- 
subtracted) saturated stellar profile from a digitised Bamberg plate (points) and 
a synthetic profile (solid line) derived from the polynomial fit to the data of King 
(1971). The synthetic profile is not a fit to the Bamberg data. The saturation 
limit of the synthetic profile was chosen so that the radius of saturation (in units 
of pixels of 5.25 arc seconds square) approximately matched the observed radius 
of saturation. 
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Figure 4: Synthetic stellar flux profiles, derived from the fit to the data of King 
(1971). Left panel - a series of profiles spanning ^4 magnitudes, plotted on 
a semi-log scale. Flux is in arbitrary units, while the horizontal axis is in arc 
seconds. Right panel - the same profiles but with a common 'saturation limit' 
imposed. The horizontal scale in this panel is now in units of 5.25 arc second 
square pixels to match the resolution of the scanned Bamberg plates. 
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Figure 5: A comparison between the observed 'radius of saturation' (open cir- 
cles) for 19 stars ranging from Tycho-2 BT magnitudes of ^13 to <~9, and the 
calculated radius of saturation (solid line) found using the King profiles. The 
flux level where saturation occurred was the only parameter adjusted to give 
the calculated curve, and was chosen on the basis of the comparison shown in 
figure 3. The observed saturation radii are based on averages from 6 'average' 
plates. The disparity at the lower radii values may in part be a consequence 
of the limited resolution of the plate scan (5.25 arc second pixels), while the 
calculated radii were derived from a 1 arc second grid and then rebinned. 



resolution of the scan, undersampling starts to reduce the observed saturation 
radius, hence Fig. 5 shows poorer agreement for the fainter stars. 
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Comparison of synthetic and aperture photometry 

The mean plate magnitudes for 19 stars, as measured on 90 plates, are 
shown as open circles in Fig. 6 plotted against their BT (Tycho-2) magnitudes. 
For comparison a 1:1 relation is shown as the solid line with upright crosses. 
The observed relationship between plate and catalogued magnitudes is reason- 
ably linear for the fainter stars, but clearly falls off (by about 1 mag) for the 
brightest stars as saturation effects become more significant. That saturation 
is the cause is clear from the aperture photometry magnitudes determined for 
the saturated, synthetic profiles shown in the right hand panel of Fig. 4, using 
a circular aperture of 20 pixels. (The same flux saturation level was used as 
that which produced the results of Fig. 5. The vertical offset is arbitrary for 
the plate and synthetic magnitudes, and has been adjusted so the data overlap 
near BT^12, which is the same zero-point as was used to plot the solid line 
in Fig. 5.) The synthetic photometry results are shown in Fig. 6 by the solid 
line with no symbols. The general trend of this line matches quite well the 
observed plate magnitudes (circles). It is relatively straightforward to derive a 
polynomial transformation from plate to standard magnitude 6 . 

Conclusion 

We have not attempted to optimise the fit to the radii of saturation (Fig. 5) 
by a more careful selection of the synthetic saturation level, or to improve the 
agreement by seeking a better match between the synthetic profile and the stellar 
profiles obtained from the Bamberg plates. The purpose of the study was to 
illustrate that the known, extended stellar radial profile, when digitised as a 
photo-positive from a photographic plate, can be used successfully for aperture 
photometry, even though (indeed, because) the stellar image is mostly saturated. 
Such a result was well known to iris photometrists (e.g. Moffat 4 , Schacfer 5 ), 
but may be less widely appreciated in the modern digital age. 
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Figure 6: Mean measured plate magnitudes for 19 stars (circles), based on aper- 
ture photometry with IRAF, from 90 digitised Bamberg plates, and calculated 
magnitudes from aperture photometry (solid line, no symbols) of the synthetic 
saturated profiles shown in the right panel of Figure 4, both plotted against 
Tycho-2 BT magnitudes. The vertical axis has an arbitrary offset (the magni- 
tude zero point) which has been adjusted to align the data near BT^12. For 
comparison a 1:1 line is also shown (solid line with + symbols.) An aperture of 
20 pixels (~ 105 arc sec) was used for both the IRAF measurements and the 
synthetic photometry. 
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